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I . INTRODUCTION 

The Penn State Finite Difference Time Domain Electromagnetic 
Scattering Code Version A is a three dimensional numerical 
electromagnetic scattering code based upon the Finite Difference 
Time Domain Technique (FDTD) . The supplied version of the code 
is one version of our current three dimensional FDTD code set. 
This manual provides a description of the code and corresponding 
results for the default scattering problem. The manual is 
organized into fourteen sections: introduction, description of 

the FDTD method, Operation, resource requirements, Version A code 
capabilities, a brief description of the default scattering 
geometry, a brief description of each subroutine, a description 
of the include file (COMMONA. FOR) , a section briefly discussing 
Radar Cross Section (RCS) computations, a section discussing the 
scattering results, a sample problem setup section, a new problem 
checklist, references and figure titles. 

II. FDTD METHOD 

The Finite Difference Time Domain (FDTD) technique models 
transient electromagnetic scattering and interactions with 
objects of arbitrary shape and/or material composition. The 
technique was first proposed by Yee [1] for isotropic, non- 
dispersive materials in 1966; and has matured within the past 
twenty years into a robust and efficient computational method. 

The present FDTD technique is capabable of transient 
electromagnetic interactions with objects of arbitrary and 
complicated geometrical shape and material composition over a 
large band of frequencies. This technique has recently been 
extended to include dispersive dielectric materials, chiral 
materials and plasmas. 

In the FDTD method, Maxwell's curl equations are discretized 
in time and space and all derivatives (temporal and spatial) are 
approximated by central differences. The electric and magnetic 
fields are interleaved in space and time and are updated in a 
second-order accurate leapfrog scheme. The computational space 
is divided into cells with the electric fields located on the 
edges and the magnetic fields on the faces (see Figure 1) . FDTD 
objects are defined by specifying dielectric and/or magnetic 
material parameters at electric and/or magnetic field locations. 

Two basic implementations of the FDTD method are widely used 
for electromagnetic analysis: total field formalism and scattered 
field formalism. In the total field formalism, the electric and 
magnetic field are updated based upon the material type present 
at each spatial location. In the scattered field formalism, the 
incident waveform is defined analytically and the scattered field 
is coupled to the incident field through the different material 
types. For the incident field, any waveform, angle of incidence 
and polarization is possible. The separation of the incident and 



scattered fields conveniently allows an absorbing boundary to be 
employed at the extremities of the discretized problem space to 
absorb the scattered fields. 
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This code is a scattered field code, and the total E and H 
fields may be found by combining the incident and scattered 
fields. Any type of field quantity (incident, scattered, or 
total) , Poynting vector or current are available anywhere within 
the computational space. These fields, incident, scattered and 
total, may be found within, on or about the interaction object 
placed in the problem space. By using a near to far field 
transformation, far fields can be determined from the near fields 
within the problem space thereby affording radiation patterns and 
RCS values. The accuracy of these calculations is typically 
within a dB of analytic solutions for dielectric and magnetic 
sphere scattering. Further improvements are expected as better 
absorbing boundary conditions are developed and incorporated. 

III. OPERATION 


Typically, a truncated Gaussian incident waveform is used to 
excite the system being modeled, however certain code versions 
also provide a smooth cosine waveform for convenience in modeling 
dispersive materials. The interaction object is defined in the 
discretized problem space with arrays at each cell location 
created by the discretization. All three dielectric material 
types for E field components within a cell can be individually 
specified by the arrays IDONE(I, J,K) , IDTWO(I, J,K) , 

IDTHRE(I, J,K) . This models arbitrary dielectric materials with m 
= /i 0 . By an obvious extension to six arrays, magnetic materials 
with m f Mo can be modeled. 


Scattering occurs when the incident wave, marched forward in 
time in small steps set by the Courant stability condition, 
reaches the interaction object. Here a scattered wave must 
appear along with the incident wave so that the Maxwell equations 
are satisfied. If the material is a perfectly conductive metal 
then only the well known boundary condition 


scat inc 

^tan = “E tan 


(1) 


must be satisfied. For a nondispersive dielectric the 
requirement is that the total field must satisfy the Maxwell 
equations in the material: 


Vx E tot 


Vx (E' nc + E scat 
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VxH tot = Vx (H inc +H scat ) 


0 E tot 
= e 

at 


+ a E 


tot 


(3) 


= e ^ (E 100 -I + g (E inc + E scat ) 

at 


(4) 


Additionally the incident wave, defined as moving unimpeded 
through a vacuum in the problem space, satisfies everywhere in 
the problem the Maxwell equations for free space 


VxE' nc 


1 dH inc 

m 0 at 


(5) 


VxH inc 


dE inc 

at 


(6) 


Subtracting the second set of equations from the first yields the 
Maxwell equations governing the scattered fields in the material: 


Vx E scat 


i aH scat 
m 0 at 


(7) 


aE 


inc 


apSMt 

,inc _ ac ‘ . _r,scat 


Vxh = (e -€ ) — — + 0 E + e 

0 at at 


+ CTE 


( 8 ) 


Outside the material this simplifies to: 


Vx e 


scat 


i da! 
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at 


(9) 


VxH scat =e 


dE 


.scat 


at 


(10) 


Magnetic materials, dispersive effects, non-linearities, 
etc., are further generalizations of the above approach. Based 
on the value of the material type, the subroutines for 
calculating scattered E and H field components branch to the 
appropriate expression for that scattered field component and 



7 


that component is advanced in time according to the selected 
algorithm. As many materials can be modeled as desired, the 
number equals the dimension selected for the flags. If materials 
with behavior different from those described above must be 
modeled, then after the appropriate algorithm is found, the 
code's branching structure allows easy incorporation of the new 
behavior. 

IV. RESOURCE REQUIREMENTS 

The number of cells the problem space is divided into times 
the six components per cell set the problem space storage 
requirements 

Storage = NC x 6 components/cell x 4 bytes /component (H) 

and the computational cost 

Operations = NC x 6 comp/cell x io ops/component x N (12) 

where N is the number of time steps desired. 

N typically is on the order of ten times the number of cells 
on one side of the problem space. More precisely for cubical 

cells it takes time steps to traverse a single cell when the 
time step is set by the Courant stability condition 

Ax 

At = Ax = cell size dimension (13) 

/7c 


The condition on N is then that 

i i 

number cells on a side 
of the problem space 


N 


10 x ( 


/7n 


C 3 ) 


NC 


(14) 


The earliest aircraft modeling using FDTD with approximately 30 
cells on a side required approximately 500 time steps. For more 
recent modeling with approximately 100 cells on a side, 2000 or 
more time steps are used. 

For (100 cell) 3 problem spaces, 24 MBytes of memory are 
required to store the fields. Problems on the order of this size 
have been run on a Silicon Graphics 4D 220 with 32 MBytes of 
memory, IBM RISC 6000, an Intel 486 based machine, and VAX 
11/785. Storage is only a problem as in the case of the 486 
where only 16 MBytes of memory was available. This limited the 
problem space size to approximately (80 cells) . 



8 


For (100 cell) problems with approximately 2000 time steps, 
there is a total of 120 x 10 9 operations to perform. The speeds 
of the previously mentioned machines are 24 MFLOPs (4 processor 
upgraded version), 10 MFLDPS, 1.5 MFLOPS, and 0.2 MFLOPs. The 
run times are then 5 3 x 10 3 seconds, 12 x 10 3 seconds, 80 x 10 3 
seconds and 600 x 10 seconds, respectively. In hours the times 
are 1.4, 3.3, 22.2 and 167 hours. Problems of this size are 
possible on all but the last machine and can in fact be performed 
on a personal computer (486) if one day turnarounds are 
permissible. 

V. VERSION A CODE CAPABILITIES 

The Penn State University FDTD Electromagnetic Scattering 
Code Version A has the following capabilities: 

1) Ability to model lossy dielectric and perfectly conducting 
scatterers . 

2) First and second order outer radiation boundary condition 
(ORBC) operating on the electric fields for dielectric or 
perfectly conducting scatterers. 

3) Near to far zone transformation capability to obtain far zone 
scattered fields. 

4) Gaussian and smooth cosine incident waveforms with arbitrary 
incidence angles. 

5) Near zone field, current or power sampling capability. 

6) Companion code for computing Radar Cross Section (RCS) . 

VI. DEFAULT SCATTERING GEOMETRY 

The code as delivered is set up to calculate the far zone 
backscatter fields for an infinitely thin, 29 cm square, 
perfectly conducting plate. The problem space size is 60 by 60 
by 49 cells in the x, y and z directions, the cells are 1 cm 
cubes, and the incident waveform is a 0-polarized Gaussian pulse 
with incidence angles of 0=45 and 0=30 degrees. The output data 
files are included as a reference along with a code (RCS3D.FOR) 
for computing the frequency domain RCS using these output data 
files. The ORBC is the second order absorbing boundary condition 
set forth by Mur [2]. 

VII. SUBROUTINE DESCRIPTION 

In the description for each subroutine, an asterisk (*) will 
be placed by the subroutine name if that particular subroutine is 
normally modified when defining a scattering problem. 
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MAIN ROUTINE 

The main routine in the program contains the calls for all 
necessary subroutines to initialize the problem space and 
scattering object (s) and for the incident waveform, far zone 
transformation, field update subroutines, outer radiation 
boundary conditions and field sampling. 

The main routine begins with the include statement and then 
appropriate data files are opened, and subroutines ZERO, BUILD 
and SETUP are called to initialize variables and/or arrays, build 
the object (s) and initialize the incident waveform and 
miscellaneous parameters, respectively. Subroutine SETFZ is 
called to intialize parameters for the near to far zone 
transformation if far zone fields are desired. 

The main loop is entered next, where all of the primary 
field computations and data saving takes place. During each time 
step cycle, the EXSFLD , EYSFLD, and EZSFLD subroutines are called 
to update the x, y, and z components of the scattered electric 
field. The six electric field outer radiation boundary 
conditions (RADE??) are called next to absorb any outgoing 
scattered fields. Time is then advanced 1/2 time step according 
to the Yee algorithm and then the HXSFLD, HYSFLD, AND HZSFLD 
subroutines are called to update the x, y, and z components of 
scattered magnetic field. Time is then advanced another 1/2 step 
and then either near zone fields are sampled and written to disk 
in DATSAV , and/or the near zone to far zone vector potentials are 
updated in SAVFZ. The parameter NZFZ (described later) in the 
common file defines the type of output fields desired. 

After execution of all time steps in the main field update 
loop, subroutine FAROUT is called if far zone fields are desired 
to compute the far zone fields and write them to disk. At this 
point, the execution is complete. 

SUBROUTINE SETFZ 

This subroutine initializes the necessary parameters 
required for far zone field computations. The code as furnished 
computes backscatter far zone fields and can compute bistatic far 
zone fields for one scattering angle (i.e. one 0 and <p angle) . 
Refer to reference [3] for a complete description of the near to 
far zone transformation. Other versions of this subroutine 
provide for multiple bistatic angles. 

SUBROUTINE SAVFZ 

This subroutine updates the near zone to far zone vector 
potentials . 
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SUBROUTINE FAROUT 

This subroutine changes the near zone to far zone vector 
potentials to far zone electric field 0 and 0 components and 
writes them to disk. 

SUBROUTINE BUILD * 

This subroutine "builds" the scattering object (s) by 
initializing the IDONE, IDTWO, and IDTHRE arrays. The 
IDONE-IDTHRE arrays are for specifying perfectly conducting and 
lossy dielectric materials. Refer to Figure 1 for a diagram of 
the basic Yee cell. For example, setting an element of the IDONE 
array at some I,J,K location is actually locating dielectric 
material at a cell edge whose center location is I+0.5,J,K. 

Thus, materials with diagonal permittivity tensors can be 
modeled. The default material type for all ID??? arrays is 0, or 
free space. By initializing these arrays to values other than 0, 
the user is defining an object by determining what material types 
are present at each spatial location. Other material types 
available for IDONE-IDTHRE are 1 for perfectly conducting objects 
and 2-9 for lossy non-magnet ic dielectrics. It is assumed 
throughout the code that all dielectric materials are non- 
magnetic (i.e. the materials have a permeability of n 0 ) . This 
subroutine also has a section that checks the ID??? arrays to 
determine if legal material types have been defined throughout 
the problem space. The actual material parameters (e and a) are 
defined in subroutine SETUP. The default geometry is a 29 cm 
square perfectly conducting plate. 

The user must be careful that his/her object created in the 
BUILD subroutine is properly formed. 

When it is important to place the object in the center of 
the problem space (to have lowest possible cross-pol scattering 
for symmetric objects) , NX etc. should be odd. This is due to 
the field locations in the Yee cell and also the placement of the 
E field absorbing boundary condition surfaces. 

If the object being modeled has curved surfaces, edges, etc. 
that are at an angle to one or more of the coordinate axes, then 
that shape must be approximately modeled by lines and faces in a 
"stair-stepped" (or stair-cased) fashion. This stair-cased 
approximation introduces errors into computations at higher 
frequencies. Intuitively, the error becomes smaller as more 
cells are used to stair-case a particular object. 

SUBROUTINE DCUBE 

This subroutine builds cubes of dielectric material by 
defining four each of IDONE, IDTWO and IDTHRE components 
corresponding to one spatial cube of dielectric material. It can 
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also be used to define thin (i.e. up to one cell thick) 
dielectric or perfectly conducting plates. Refer to comments 
within DCUBE for a description of the arguments and usage of the 
subroutine. 

SUBROUTINE SETUP * 

This subroutine initializes many of the constants required 
for incident field definition, field update equations, outer 
radiation boundary conditions and material parameters. The 
material parameters e and a are defined for each material type 
using the material arrays EPS and SIGMA respectively. The array 
EPS is used for the total permittivity and SIGMA is used for the 
electric conductivity. These arrays are initialized in SETUP to 
free space material parameters for all material types and then 
the user is required to modify these arrays for his/her 
scattering materials. Thus, for the lossy dielectric material 
type 2, the user must define EPS(2) and SIGMA(2) . The remainder 
of the subroutine computes constants used in field update 
equations and boundary conditions and writes the diagnostics 
file. 

SUBROUTINE EXSFLD 

This subroutine updates all x components of scattered 
electric field at each time step except those on the outer 
boundaries of the problem space. IF statements based upon the 
IDONE array are used to determine the type of material present 
and the corresponding update equation to be used. These 
scattered field equations are based upon the development given in 
[4]. 

SUBROUTINE EYSFLD 

This subroutine updates all y components of scattered 
electric field at each time step except those on the outer 
boundaries of the problem space. IF statements based upon the 
IDTWO array are used to determine the type of material present 
and the corresponding update equation to be used. 

SUBROUTINE EZSFLD 

This subroutine updates all z components of scattered 
electric field at each time step except those on the outer 
boundaries of the problem space. IF statements based upon the 
IDTHRE array are used to determine the type of material present 
and the corresponding update equation to be used. 

SUBROUTINES RADEYX , RADEZX, RADEZY, RADEXY, RADEXZ and RADEYZ 

These subroutines apply the outer radiation boundary 
conditions to the scattered electric field on the outer 



boundaries of the problem space. 

SUBROUTINE HXSFLD 

This subroutine updates all x components of scattered 
magnetic field at each time step. The standard non-magnetic 
update equation is used. 

SUBROUTINE HYSFLD 

This subroutine updates all y components of scattered 
magnetic field at each time step. The standard non-magnetic 
update equation is used. 

SUBROUTINE HZSFLD 

This subroutine updates all z components of scattered 
magnetic field at each time step. The standard non-magnetic 
update equation is used. 

SUBROUTINE DATSAV * 

This subroutine samples near zone scattered field quantities 
and saves them to disk. This subroutine is where the quantities 
to be sampled and their spatial locations are to be specified and 
is onl y c a lle d if near zone f ields only are de sired or if both 
near, an d far zone fields are desired . Total field quantities can 
also be sampled. See comments within the subroutine for 
specifying sampled scattered and/or total field quantities. When 
sampling magnetic fields, remember the <St/2 time difference 
between E and H when writing the fields to disk. Sections of 
code within this subroutine determine if the sampled quantities 
and the spatial locations have been properly defined. 

FUNCTIONS EXI, EYI and EZI 

These functions are called to compute the x, y and z 
components of incident electric field. The functional form of 
the incident field is contained in a separate function SOURCE. 

FUNCTION SOURCE * 

This function contains the functional form of the incident 
field. The code as furnished uses the Gaussian form of the 
incident field. An incident smooth cosine pulse is also 
available by uncommenting the required lines and commenting out 
the Gaussian pulse. Thus, this function need only be modified if 
the user changes the incident pulse from Gaussian to smooth 
cosine. A slight improvement in computing speed and 
vectorization may be achieved by moving this function inside each 
of the incident field functions EXI, EYI and so on. 
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FUNCTIONS DEXI, DEYI and DEZI 

These functions are called to compute the x, y and z 
components of the time derivative of incident electric field. 

The functional form of the incident field is contained in a 
separate function DSRCE . 

FUNCTION DSRCE * 

This function contains the functional form of the time 
derivative of the incident field. The code as furnished uses the 
time derivative of the Gaussian form of the incident field. A 
smooth cosine pulse time derivative is also available by 
uncommenting the required lines and commenting out the Gaussian 
pulse. Thus, the function need only be modified if the user 
changes from the Gaussian to smooth cosine pulse. Again, a 
slight improvement in computing speed and vectorization may be 
achieved by moving this function inside each of the time 
derivative incident field functions DEXI, DEYI and so on. 

SUBROUTINE ZERO 

This subroutine initializes various arrays and variables to 

zero. 

VIII. INCLUDE FILE DESCRIPTION (COMMONA. FOR) * 

The include file, COMMONA. FOR, contains all of the arrays 
and variables that are shared among the different subroutines. 
This file will require the most modifications when defining 
scattering problems. A description of the parameters that are 
normally modified follows. 

The parameters NX, NY and NZ specify the size of the problem 
space in cells in the x, y and z directions respectively. For 
problems where it is crucial to center the object within the 
problem space, then NX, NY and NZ should be odd. The parameter 
NTEST defines the number of near zone quantities to be sampled 
and NZFZ defines the field output format. Set NZFZ=0 for near 
zone fields only, NZFZ=1 for far zone fields only and NZFZ=2 for 
both near and far zone fields. Parameter NSTOP defines the 
maximum number of time steps. DELX, DELY, and DELZ (in meters) 
define the cell size in the x, y and z directions respectively. 
The 9 and <p incidence angles (in degrees) are defined by THINC 
and PHINC respectively and the polarization is defined by ETHINC 
and EPHINC. ETHINC=1 . 0 , EPHINC=0.0 for ©-polarized incident 
field and ETHINC=0.0, EPHINC=1.0 for 0-polarized incident fields. 
Parameters AMP and BETA define the maximum amplitude and the e 
temporal width of the incident pulse respectively. BETA 
automatically adjusts when the cell size is changed and normally 
should not be changed by the user. The far zone scattering 
angles are defined by THETFZ and PHIFZ. The code as furnished 
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perforins backscatter computations, but these parameters could be 
modified for a bistatic computation. 

IX. RCS COMPUTATIONS 

A companion code, RCS3D.FOR, has been included to compute 
RCS versus frequency. It uses the file name of the FDTD far zone 
output data (FZ0UT3D.DAT) and writes a data file of far zone 
electric fields versus time (FZTIME.DAT) and RCS versus frequency 
(3DRCS.DAT) . The RCS computations are performed up to the 10 
cell/l 0 frequency limit. Refer to comments within this code for 
further details. 

X. RESULTS 

As previously mentioned, the code as furnished models an 
infinitely thin, 29 cm square, perfectly conducting plate and 
computes backscatter far zone scattered fields at angles of 0=45 
and 0=30 degrees. 

Figures 2-3 shows the co-polarized far zone electric field 
versus time and the co-polarized RCS for the 29 cm square 
perfectly conducting plate. 

Figures 4-5 shows the cross-polarized far zone electric 
field versus time and the cross-polarized RCS for the 29 cm 
square perfectly conducting plate. 

XI. SAMPLE PROBLEM SETUP 

The code as furnished models an infinitely thin, 29 cm 
square, perfectly conducting plate and computes backscatter far 
zone scattered fields at angles of 0=45 and 0=30 degrees. The 
corresponding output data files are also provided, along with a 
code to compute RCS using these data files. In order to change 
the code to a new problem, many different parameters need to be 
modified. A sample problem setup will now be discussed. 

Suppose that the problem to be studied is RCS backscatter 
versus frequency from a 28 cm by 31 cm perfectly conducting plate 
with a 3 cm dielectric coating with a dielectric constant of 4e 0 
using a 0-polarized field. The backscatter angles are 0=30.0 and 
0=60.0 degrees and the frequency range is up to 3 Ghz. 

Since the frequency range is up to 3 Ghz, the cell size must 
be chosen appropriately to resolve the field IN ANY MATERIAL at 
the highest frequency of interest. A general rule is that the 
cell size should be 1/10 of the wavelength at the highest 
frequency of interest. For difficult geometries, 1/20 of a 
wavelength may be necessary. The free space wavelength at 3 GHz 
is 1 0 =10 cm and the wavelength in the dielectric coating at 3 GHz 
is 5 cm. The cell size is chosen as 1 cm, which provides a 
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resolution of 5 cells/A in the dielectric coating and 10 cells/A 0 
in free space. Numerical studies have shown that choosing the 

cell size < 1/4 of the shortest wavelength in any material is 
the practical lower limit. Thus the cell size of 1 cm is barely 
adequate. The cell size in the x, y and z directions is set in 
the common file through variables DELX, DELY and DELZ. Next the 
problem space size must be large enough to accomodate the 
scattering object, plus at least a five cell boundary (10 cells 
is more appropriate) on every side of the object to allow for the 
far zone field integration surface. It is advisable for plate 
scattering to have the plate centered in the x and y directions 
of the problem space in order to reduce the cross-polarized 
backscatter and to position the plate low in the z direction to 
allow strong specular reflections multiple encounters with the 
ORBC. A 10 cell border is chosen, and the problem space size is 
chosen as 49 by 52 by 49 cells in the x, y and z directions 
respectively. As an initial estimate, allow 2048 time steps so 
that energy trapped within the dielectric layer will radiate. 

Thus parameters NX, NY and NZ in COMMONA.FOR would be changed to 
reflect the new problem space size, and parameter NSTOP is 
changed to 2048. If all transients have not been dissipated 
after 2048 time steps, then NSTOP will have to be increased. 
Truncating the time record before all transients have dissipated 
will corrupt frequency domain results. Parameter NZFZ must be 
equal to l since we are interested in far zone fields only. To 
build the object, the following lines are inserted into the BUILD 
subroutine: 

C 

C BUILD THE DIELECTRIC SLAB FIRST 

C 

ISTART=11 
JSTART=11 
KSTART=1 1 
NXWIDE=2 8 
NYWIDE=3 1 
NZWIDE=3 
MTYPE=2 

CALL DCUBE ( ISTART , JSTART , K START , NXWIDE , NYWIDE , NZWIDE , MTYPE) 
C 

C BUILD PEC PLATE NEXT 

C 

ISTART=11 

JSTART=11 

KSTART=11 

NXWIDE=2 8 

NYWIDE=3 1 

NZWIDE=0 

MTYPE=1 

CALL DCUBE ( ISTART , JSTART , KSTART , NXWIDE , NYWIDE , NZWIDE , MTYPE) 
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The PEC plate is built last on the bottom of the dielectric 
slab to avoid any air gaps between the dielectric material and 
the PEC plate. In the common file, the incidence angles THINC 
and PHINC have to be changed to 30.0 and 60.0 respectively, the 
cell sizes ( DELX , DELY, DELZ) are set to 0.01, and the 
polarization is set to ETHINC=1.0 and EPHINC=0.0 for 0-polarized 
fields. Since dielectric material 2 is being used for the 
dielectric coating, the constitutive parameters EPS (2) and 
SIGMA(2) are set to 4e 0 and 0.0 respectively, in subroutine 
SETUP. This completes the code modifications for the sample 
problem. 

XII. NEW PROBLEM CHECKLIST 

This checklist provides a quick reference to determine if 
all parameters have been defined properly for a given scattering 
problem. A reminder when defining quantities within the code: 
use MKS units and specify all angles in degrees. 

COMMONA . FOR : 

1) Is the problem space sized correctly? (NX, NY, NZ) 

2) For near zone fields, is the number of sample points correct? 
(NTEST) 

3) Is parameter NZFZ defined correctly for desired field 
outputs? 

4) Is the number of time steps correct? (NSTOP) 

5) Are the cell dimensions (DELX, DELY, DELZ) defined correctly? 

6) Are the incidence angles (THINC, PHINC) defined correctly? 

7) Is the polarization of the incident wave defined correctly 
(ETHINC, EPHINC)? 

8) For other than backscatter far zone field computations, are 
the scattering angles set correctly? (THETFZ, PHIFZ) 

SUBROUTINE BUILD: 

1) Is the object completely and correctly specified? 

SUBROUTINE SETUP: 

1) Are the constitutive parameters for each material specified 
correctly? (EPS and SIGMA) 
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FUNCTIONS SOURCE and DSRCE: 

1) If the Gaussian pulse is not desired, is it commented out and 

the smooth cosine pulse uncommented? 

SUBROUTINE DATSAV : 

1) For near zone fields, are the sampled field types and spatial 

locations correct for each sampling point? (NTYPE, IOBS, JOBS, 

KOBS) 
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XIV. FIGURE TITLES 

Fig. 1 Standard three dimensional Yee cell showing placement 
of electric and magnetic fields. 

Fig. 2 Co-polarized far zone scattered field versus time for 
29 cm square perfectly conducting plate. 

Fig. 3 Co-polarized RCS versus frequency for 29 cm square 
perfectly conducting plate. 

Fig. 4 Cross-polarized far zone scattered field versus time 
for 29 cm square perfectly conducting plate. 

Fig. 5 Cross-polarized RCS versus frequency for 29 cm square 
perfectly conducting plate. 
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